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[0003] 

Background of Invention 
The Field of the Invention 

[0004] Thj S invention relates to the field of robust control engi- 
neering, specifically to the design of robust controllers 
that can make the overall system performance less sensi- 



tive to the model imperfection of system under control. 

Some of the related fields include 700/28, 700/29, and 

700/38. 
Background Information 

[0005] Robustness is an important issue in many engineering 

disciplines. In the field of control engineering, robustness 
against disturbances and model uncertainty is at the heart 
of control practice. The key part of robust feedback con- 
trol has been focusing on the effects of plant modeling 
uncertainty on stability. This is in contrast with optimal 
control, which usually deals with disturbance rejection 
while assuming the absence of modeling error. Robust 
optimal control is a combination of the two, which can 
keep certain level of disturbance rejection despite the 
presence of modeling error. If a control design is opti- 
mized only for a nominal model without considering pos- 
sible modeling error, it can perform far from expectation, 
since the actual system dynamics is not represented ex- 
actly by the nominal model. Conversely, it a control de- 
sign is optimized only for maximizing a stability margin 
without considering the magnitude and inherent structure 
of the actual modeling error, the resulting system perfor- 
mance can be over conservative, since the actual system 



dynamics variation is not as large as modeled in the con- 
trol design stage. In addition, if performance specifica- 
tions can not be directly specified in the control design 
problem formulation, it is impractical to expect the result- 
ing system performance can be meet the specifications. 
[0006] For designing a controller satisfying some robust stability 
measures, the D-K iteration method is a widely-adopted 
approach. D-K iteration was first proposed in the early 
80"s. One of the key publications is [J. Doyle, "Analysis of 
feedback systems with structured uncertainties," IEE Pro- 
ceedings, 129(6), part D, Nov. 1982]. Some prior-art pub- 
lications that explain D-K iteration in details of includes 
[S. Skogestad et al, Multivariable Feedback Control, John 
1996][K. Zhou, Essentials of Robust Control, Prentice Hall, 
1998][G. Balas et al, u-Analysis and Synthesis Toolbox, 
version 3, The Mathworks, Inc., 1998][R. Sanchez-Pena et 
al, Robust Systems Theory and Applications, 1998]. It 
starts with defining a general control problem formulation 
as shown in Fig. I, where P is a generalized plant model 
containing a nominal plant model, one or multiple model 
uncertainty weighting filters to describe the variation of 
plant uncertainties at different frequencies, and one or 
multiple selected performance weighting filters to adjust 



the performance level of the system; K is a controller to be 
designed; and A represents the dynamics of plant varia- 
tion. For some system performance measures, such as the 
H norm, the robust optimal control problem can be for- 

00 

mulated as a augmented robust stability as shown in Fig 

IB, where A represents some fictitious uncertainty for the 
p 

performance channel, from the exogenous input w to the 
exogenous output z. System no, which is same with sys- 
tem 120 denoted by A , is a resulting uncertainty system 

a 

by augmenting A to A. In D-K iteration, a frequency- vary 

p 

block-diagonal scaling D in Fig. 1C is used in its D-step to 
exploit the structure of the augmented uncertainty A^. 
Thus it can result in a robustly stable system without 
trad ing-off too much performance. 
[0007] a typical D-K iteration procedure can be shown with the 
example flowchart in Fig. 3. In step 310, a generalized 
plant model is defined, which represent the interconnec- 
tions between a provided nominal plant model, the mod- 
els of one or multiple provides performance weighting fil- 
ters, and the structure of plant uncertainties. In step 320, 
a parametric model of the uncertainty weighting filter de- 
noted by W is provided. It can be an assumed variation, 

del 

or it can be estimated from a set of experimental data 



which captures the actual variation of the system. After 
step 310 and step 320, a general control problem is for- 
mulated with all modeling data required for the following 
steps. Step 330 initializes the iteration by assigning a ini- 
tial value D q the block-diagonal uncertainty scaling D in 
Fig. 1C. A common choice of D q is a block-diagonal iden- 
tity matrix. In step 330 and step 350, a robustness mea- 
sure, usually an upper bound of a structured singular val- 
ued denoted by u, is minimized with respect to the con- 
troller K, while holding D fixed. Explanations of u and its 
upper bound related to H norm are explained in many of 

00 

the prior art publications, such as [A. Packard et al "Lin- 
ear, multivariable robust control with a u perspective," 
ASME J. of Dynamics, measurement, and Control, vol. 115 
Jun. 1993, pp. 426-438]. A robustness measure can also 
be specified with L norm other than u, as described in [A. 
Dahleh et al, Control of Uncertain Systems, 1996]. In com- 
mon practices, u is selected as the robustness measure. 
With this selection, in step 330 and step 3 5 0, K is solved 
as an optimal H controller, commonly by using state- 

00 

space approaches such as J. Doyle et al "State-space solu- 
tions to standard H and H control problems," IEEE 
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Transactions on Automatic Control, 34(8): 831-847, 



1989][P. Gahinet et al "A linear matrix inequality approach 
to H control," International J. Robust Nonlinear Control, 

00 

(4)421-448, 1994]. The first is a state-space approach 
based on a set of Riccati equations, and the second is a 
generalized of the first by linear matrix equalities. In step 
340, the D is commonly computed by first optimizing the 
frequency response of D at each of a set of selected fre- 
quencies independently, based on frequency gridding op- 
timization. However, the state-space approaches in step 
350 require parametric models as input data, therefore in 
step 3 4 0, the frequency response of D needs to be curve- 
fitted to be used in step 350. Similarly, in step 320 para- 
metric models of the plant uncertainty weight W de| needs 
to be provided by curve-fitting the experimental data, 
usually the frequency response of W . Step 340, common 

del 

referred as D-step, and step 350, common referred as K- 
step, are iterated until a stopping criterion is met. Com- 
mon stopping criterions includes the selected robustness 
measure being smaller than a desired value, or the ro- 
bustness measure not decreasing with more iterations. 
[0008] There are several issues with the D-K iteration method. 
First, the requirement of curve-fitting in the step 320 for 
Wdel and step 340 for D can be problematic. When the fre- 



quency response of Wdel or D is complicated, it is hard to 
curve-fit an accurate parametric model. This can easily 
cause the designed controller be over conservative. In ad- 
dition, the error of the D curve-fitting in each D-K itera- 
tion can cause the upper bound of u be increasing in 
practice. It is also possible that when attempting to use a 
higher-order model to improve the curve-fitting fidelity, 
numerical problems such as ill-conditioning can occur. 
[0009] Second, it is hard to directly incorporate multiple, espe- 
cially time domain specifications into the D-K iteration 
design method. Time domain specifications, such as the 
rise time, the peak step response overshoot, and the peak 
control effort, etc., are very important for practical control 
systems. However, D-K iteration, along with most other 
well-known control design approaches (e.g. P-l-D control, 
lead-lag control, and LQG method,) can not take these 
specifications into account directly. The common practice 
is trying to tune "weighting filters" to meet specifications 
indirectly. It can be a very long and iterative process if 
there are many specifications imposed together. In addi- 
tion, the use of a single Hoo norm to measure the system 
performance is usually over-simplified. Practical perfor- 
mance specifications are usually specified on different in- 



put-output channels, thus using a single system norm to 
approximate them can easily lead to performance over- 
conservatism. The state-space approach based on linear 
matrix equalities does allow multiple controller design 
constrains to be incorporated together, however in many 
cases some constraints on the decision variables need to 
be imposed such that the overall optimization problem is 
still convex, as shown in [P. Gahinet, et al. LMI Control 
Toolbox. The Mathworks, Inc, 1995]. Therefore it can eas- 
ily lead to performance over-conservatism. 
[0010] since late 80"s and early 90"s, Q-parameterization design 
method has been widely applied to design controllers that 
can incorporate multiple time-domain and frequency-do- 
main performance specifications. Some of the prior-art 
publications include [B. Rafaely et al, "H /H active con- 
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trol of sound in a headrest: design and implementation," 
IEEE Trans. Control System Technology, vol. 7, no. l,Jan. 
1999][P. Titterton, "Practical method for constrained-op- 
timization controller design: H or H optimization with 
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multiple H and/or H constraints," IEEE Proceedings of 
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ASILO 1996][P. Titterton, "Practical multi-constraint H 

00 

controller synthesis from time-domain data," International 
J. of Robust and Nonlinear Control, vol. 6, 413-430, 



1996][S. P. Wu et al, "FIR filter design via spectral factor- 
ization and convex optimization," in Applied Computa- 
tional Control, Signal and Communications, 1997][K. Tsai 
et al, "DQIT: u-synthesis without D-Scale Fitting," Ameri- 
can Control Conference 2002, pp. 493-498][S. Boyd et al, 
Linear Controller Design: Limits of Performance. Prentice- 
Hall, 1991] and [S. Boyd et al, "A new CAD method and 
associated architectures for linear controllers," IEEE Trans- 
actions on Automatic Control, vol.33, p. 268, 1988]. Fig. 
ID and Fig. IE show that after Q-parameterization, the 
new free controller design parameter becomes Q. The 
method starts with transforming the generalized control 
design problem formulation in P and K of Fig 1A to a new 
form as shown in Fig. 1 £, in general using a technique 
known as Q-parameterization or Youla-parameterization. 
Fig. ID shows the structure of Q-parameterization, where 
a stabilizer J stabilize the system, and the equalizer Q is 
used to adjust the system response without causing insta- 
bility, as long as Q itself is stable. The system 210, which 
is the combination of P and J, is equivalent to the system 
N in Fig. 1 E. Once the design problem has been trans- 
formed to the form in Fig. 1 E, it can be shown that in the 
frequency domain, the exogenous output z is related to 



the exogenous input w as z = (N +N QN )w, where N 

zw zy uw zw 

is the sub-part of N transferring from w to z, is the 
sub-part of N transferring from y to z, and N is the 

u w 

sub-part of N transferring from w to u. The important as- 
pect is that at each frequency, the equalized, closed-loop 
transfer matrix from w to z is (N +N QN ), which is 

zw zy uw 

convex in terms of the frequency response of Q at the 
same frequency, the equalizer to be designed. Therefore, 
frequency-shaping specifications and the tradeoff be- 
tween different input-output channels can be specified as 
convex objectives and convex constraints when formulat- 
ing a multi-objective optimization problem. In addition, 
the frequency response data of N , N , and N can be 

zw zy uw 

incorporated without curve-fitting them first. However, 
plant uncertainties are commonly approximated with a H 

00 

norm constraint without exploiting their inherent conser- 
vatism. This can easily lead to performance over- 
conservatism. 

1 ] Therefore, what is desired is combine the strength of the 
D-K iteration method which is the capability to exploit the 
structure of uncertainties, and the strength of the Q- 
parameterization method which is the capability to incor- 
porate multiple performance objectives. The invention in- 



volves with applying Q-parameterization to Fig. 1A, and 
thus transforming the robust control design problem for- 
mulation from the A-P-K form in Fig 1A, to the A-N-Q 
form in Fig. 2B. The embodiments of the invention provide 
the capability to synthesize a controller with multiple per- 
formance objectives, while satisfying a robustness mea- 
sure which considers the inherent structure of uncertain- 
ties. In addition, the control optimization problem can be 
formulated based on frequency gridding, such that the 
numerical problems and conservatisms associated to the 
curve-fitting of uncertainty scaling in step 340 and the 
curve-fitting of plant uncertainty weight in step 320 can 
be avoided. 

[0012] Some of the related prior-work close to this invention in- 
clude [A. Lanzon et al "A Frequency Domain Optimisation 
Algorithm for Simultaneous Design of Performance 
Weights and Controllers in mu-Synthesis", Proceedings of 
the 38th IEEE Conference on Decision and Control, Vol. 5, 
pp. 4523-4528, Phoenix, AZ, USA, Dec 1999][A. Lanzon, 
"A State-Space Algorithm for the Simultaneous Optimisa- 
tion of Performance Weights and Controllers in mu- 
Synthesis", Proceedings of the 39th IEEE Conference on 
Decision and Control, Vol. 1, pp. 611-616, Sydney, Aus- 



tralia, Dec 2000 ][A Lanzon, Ph.D. Thesis: "Weight Selec- 
tion in Robust Control: An Optimisation Approach", Uni- 
versity of Cambridge, UK, Oct 2000]. In short these prior- 
art publications are denoted by [Lanzon CDC1999], 
[Lanzon CDC2000], [Lanzon PhD 2000]. In [Lanzon 
CDC1999], a control design problem is formulated equiv- 
alently to the A -N-Q form in Fig. 2C to solve a u- 

a 

synthesis problem by iterating a D-step without D-fitting, 
and a Q-parameterization design step to minimize a ro- 
bustness measure. Both the D-step and the Q-step are 
formulated based on frequency gridding optimization. 
However, it does not incorporate multiple performance 
design objectives into the iteration. The potential advan- 
tage of performing the iteration without D-fitting is not 
addressed. There is only one sentence in [Lanzon PhD 
2000, page 59] "Performance weights and D-scales are 
found and used pointwise in frequency and hence need 
not be fitted with stable minimum-phase transfer function 
matrices", without actually indicating or showing the ben- 
efits of performing the iteration without curve-fitting, as 
opposed to the publications claimed in the invention [K. 
Tsai and H. Hindi, "DQIT: u-synthesis without D-Scale Fit- 
ting," American Control Conference 2002, pp. 493-498] 



and [K. Tsai, Design of Feedforward and Feedback Con- 
trollers by Signal Processing and Convex Optimization 
Techniques, chapter 2, chapter 3, and page 129-130.] In 
fact, in [Lanzon CDC1999] it is implied by the author that 
the frequency-gridding optimization is not the preferred 
approach. Therefore in [Lanzon CDC2000], a state-space 
approach is proposed to perform the synthesis without 
frequency-gridding, and notably the A-N-Q formulation is 
abandoned, and the A-P-K formulation is used. Although 
it claims the new state-space approach can incorporate 
"other closed-loop objectives such as regional pole place- 
ment, H 2 norm minimization, etc", it is well-known that 
this is equivalent with the previously-mentioned state- 
space Hoo control solution using linear matrix inequities, 
which can easily lead to conservatism where additional 
constraints on its decision variables are imposed in order 
to preserve the convexity of the control optimization 

problem. 
Summary of Invention 

[0013] a method for design of a multi-objective least conserva- 
tive robust controller to control a plant or a process which 
may be modeled imperfectly is disclosed. It comprises a 
robust analysis step and a robust multi-objective con- 



troller synthesis step using Q-parameterization control 
design technique. The main advantages of this method 
include: 1. It allows for the tradeoff between multiple 
time-domain and frequency-domain performance objec- 
tives while keeping a robustness measure under a least 
conservative level; 2. It does not require the designer to 
provide curve-fitted parametric models for the uncertainty 
scaling and the uncertainty weight, therefore the potential 
numerical problems and performance conservatism due to 
curve-fitting can be avoided. 
[0014] | n one embodiment of the invention, the K-step of stan- 
dard D-K iteration for mu-synthesis is replaced by a Q- 
parameterization control design step. The Q-step opti- 
mization problem formulation comprises a standard ro- 
bustness measure and one or a plurality of other perfor- 
mance measures. During the iteration, the Q-step opti- 
mization problem formulation can be changed. In another 
embodiment, a controller satisfying a level of robustness 
measure is first found. Then, a Q-parameterization con- 
trol design step is performed, such that one or plurality of 
the other performance measures are optimized, while still 
satisfying a level of robustness measure which is the same 
with, or slightly traded-off from the previous level of ro- 



bustness measure. In all embodiments of the invention, if 
the robustness measure in the Q-step is formulated based 
on frequency-gridding, the problematic D-step curve fit- 
ting process in standard D-K iteration can be avoided. In 
addition, a least-conservative non-parametric plant un- 
certainty weights can incorporated directly without curve- 
fitting. Therefore the difficulties of curve-fitting and the 
conservativeness due to curve-fitting in standard D-K it- 
eration can both be eliminated. 
[0015] Although many details have been included in the descrip- 
tion and the figures, the invention is defined by the scope 
of the claims. Only limitations found in those claims apply 

to the invention. 
Brief Description of Drawings 

[0016] Figure 1A is a prior-art block-diagram showing the gen- 
eral control problem formulation for a system with model- 
ing uncertainties. 

[0017] Figure I £ is a prior-art block-diagram showing the per- 
formance uncertainties of a system can be augmented to 
its modeling uncertainties in the linear robust control 
framework 

[0018] Figurel C is a prior-art block-diagram showing the use of 
a block-diagonal scaling D on the A-P-K control problem 



formulation, to reduce the conservatism of a robustness 
measure by exploiting the structure of the uncertainty, as 
in the prior-art D-K iteration method. 

[0019] Figure ID is a prior-art block diagram showing a general- 
ized control design problem with Q-parameterization. 

[0020] Figure I £ is a prior-art block showing an equivalent gen- 
eralized control design problem with Fig ID, with Q- 
parameterization. 

[0021] Figure 2 A is a block-diagram showing the application of 
Q-parameterization on the generalized plant model P. 

[0022] Figure 2B is a block-diagram showing after Q- 

parameterization, the new generalized plant model is N, 
and the new control design parameter is Q. 

[0023] Figure 2C is a block-diagram for showing the use of a 

block-diagonal scaling D on the A-N-Q control problem 
formulation to reduce the conservatism by exploiting the 
structure of the uncertainty of a robustness measure 

[0024] Figure 3 is a prior-art flowchart of the standard D-K itera- 
tion method. 

[0025] Figure 4 is a control design flowchart as one example of 
one embodiment of the invention, which enables synthe- 
sizing multiple performance objectives by replacing the K 
step in Fig. 3 of D-K iteration with a Q-step optimization 



problem formulation, and problem and conservatism due 
to the curve-fitting of the D-scaling and plant uncertainty 
weight of D-K iteration can be avoided by formulating the 
optimization problem based on frequency-gridding. 

[0026] Figure 5 is a control design flowchart as one example of 
another embodiment of the invention, where a controller 
satisfying a level of robustness measure is first found, fol- 
lowed by a Q-parameterization control design step which 
optimizes one or multiple design objectives while still 
closely satisfying the level of robustness measure. 

[0027] Figure 6 is a control design flowchart as one example of 

applying the embodiment shown in Fig. 5 to compare with 

a prior-art multi-objective control design method without 

considering the structure of the uncertainty. 
Detailed Description 

An Overview of an Embodiment of the Invention 

[0028] a method for fast design of frequency-shaping multi- 
objective equalizers is disclosed. It comprises a robust 
analysis step and a robust multi-objective controller syn- 
thesis step using Q-parameterization control design tech- 
nique. In all embodiments of the invention, if the robust- 
ness measure in the Q-step is formulated based on fre- 



quency-gridding, the problematic D-step curve fitting 
process in standard D-K iteration can be avoided. In addi- 
tion, a least-conservative non-parametric plant uncer- 
tainty weights can incorporated directly without curve- 
fitting. Therefore the difficulties of curve-fitting and the 
conservativeness due to curve-fitting in standard D-K it- 
eration can both be eliminated. 

[0029] | n one embodiment of the invention, the K-step of stan- 
dard D-K iteration for mu-synthesis is replaced by a Q- 
parameterization control design step. The Q-step opti- 
mization problem formulation comprises a standard ro- 
bustness measure and one or a plurality of other perfor- 
mance measures. During the iteration, the Q-step opti- 
mization problem formulation can be changed. In another 
embodiment, a controller satisfying a level of robustness 
measure is first found. Then, a Q-parameterization con- 
trol design step is performed, such that one or plurality of 
the other performance measures are optimized, while still 
satisfying a level of robustness measure which is the same 
with, or slightly traded-off from the previous level of ro- 
bustness measure. 

[0030] An example design flowchart in Fig. 4 is provided for the 
first embodiment. Another example design flowchart in 



Fig. 5 is provided for the second embodiment. Comparing 
with the prior-art design flow chart in Fig. 3, the design 
steps that appear in both the prior-art and the invention 
will only be briefly explained since their details can be 
found in several prior-art references. 
[0031] a practical design example following the example design 
flowchart in Fig. 6 is described to demonstrate the effec- 
tiveness of the algorithms disclosed in the invention. 
Embedding Nominal Performance Specifications in #-synthesis 

[0032] pig 4. shows an example design flowchart for one embod- 
iment of the invention, where the control optimization 
step 450 incorporate standard robustness measure and 
one or multiple of performance objectives. The method 
starts with step 410 to define a generalized plant model, 
as in the prior-art step 310. Then Q-parameterization is 
performed in step 430, such that the problem formulation 
is in the form of Fig. 2B and Fig 2C. Then in step 440, the 
frequency response of the block-diagonal uncertainty 
scaling D is optimized frequency-by-frequency, as in the 
prior-art step 340. If in step 450, the designer chooses the 
robustness measure to be formulated based on a fre- 
quency-by-frequency gridding optimization formulation, 
then in step 440, no curve-fitting to the frequency re- 



sponse of D is required; and in step 420, the frequency re- 
sponse data of the uncertainty weights can be provided 
without curve-fitting it to a parametric model. Some 
prior-art publication explaining the procedure to perform 
this frequency gridding optimization can be found in [B. 
Rafaely et al, "H /H active control of sound in a headrest: 
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design and implementation," IEEE Trans. Control System 
Technology, vol. 7, no. 1, Jan. 1999][P. Titterton, "Practical 
method for constrained-optimization controller design: 
or H optimization with multiple H and/or H con- 

oo 2 oo 

straints," IEEE Proceedings of ASILO 1996]][A. Lanzon et al 
"A Frequency Domain Optimisation Algorithm for Simulta- 
neous Design of Performance Weights and Controllers in 
mu-Synthesis", Proceedings of the 38th IEEE Conference 
on Decision and Control, Vol. 5, pp. 4523-4528, Phoenix, 
AZ, USA, Dec 1999] [K. Tsai et al, "DQIT: u-synthesis 
without D-Scale Fitting," American Control Conference 
2002, pp. 493-498]. If the robustness measure is formu- 
lated based on non-frequency gridding approaches, then 
parametric models in step 420 and step 440 are still re- 
quired. In step 450, the control optimization formulation 
no only include the robustness measure, but also one or 
multiple performance objectives. The trade-off between 



these performance objectives and the improvement of the 
robustness measure can be adjusted by modified the con- 
trol optimization formulation in each iteration. As an ex- 
ample, a control optimization problem can be formulated 
as: minimize {(l-rho)*sigma(DH(Q)D _1 ))+rho*f o (H(Q)) } for 
a set of selected frequencies, with respective to the free 
controller design parameter Q, subject to: {f (H(Q))<0} 

k 

where k is a nonnegative integer. Here H(Q) represents 
system 130, sigma(DH(Q)D -1 )) represents the upper bound 
of u as the robustness measure, f is a performance ob- 
jective such as maximum control effort, f represents one 

k 

or multiple performance constraints, such as the noise 
amplification of one of the input-output channels, and rho 
is a weighting factor between 0 and 1, which can be ad- 
justed at each iteration to enforce the optimization 
weights more on the robustness measure or the perfor- 
mance objective f . Many other variations from this exam- 
ple optimization formulation are possible. The D-step 440 
and Q-step 450 are iterated until a decision is made to 
stop the iteration, commonly when as performance re- 
quirement has been met, or when there is no more per- 
formance improvement with more iterations. 
Reducing Robust Stability Conservatism of Multi-objective 
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[0033] pig 5. shows an example design flowchart for another em- 
bodiment of the invention, where step 550 shows one or 
multiple performance objectives can be simultaneously 
optimized while satisfying a predetermined robustness 
measure. The method starts with step 5 10 to define a 
generalized plant model, as in the prior-art step 310. Then 
Q-parameterization is performed in step 5 30, such that 
the problem formulation is in the form of Fig. 2B and Fig 
2C. In step 540 a Q-step and a D-step iteration is per- 
formed until a robustness measure is met. In fact, in one 
variation step 530 and step 540 can be replaced by the 
prior-art D-K iteration, followed by a Q-parameterization 
step before step 550. In step 550, the control optimization 
formulation includes one or multiple performance objec- 
tives or constraints, and the robustness measure intend- 
ing to keep the level of robustness measure obtained by 
step 540. As an example, suppose from step 540 an upper 
bound of u is found to be gammaO. A control optimization 
problem can be formulated as: minimize {f Q (Q)for a set of 
selected frequencies, with respective to the free controller 
design parameter Q, subject to {sigma(DH(Q)D 
_1 ))<gammal, and f (Q)<0, k=l,2,3...}. Here H(Q) repre- 



sents system 130, sigma(DH(Q)D -1 )) represents the upper 
bound of u as the robustness measure, , f is a perfor- 
mance objective such as maximum control effort, f rep- 

k 

resents one or multiple performance constraints, such as 
the noise amplification of one of the input-output chan- 
nels, and gammal is the same or slightly adjusted to be 
larger than gammaO obtained from step 540, such that the 
optimization process can have enough feasible set in or- 
der to satisfy the performance constraints f (Q)<0. Many 

k 

other variations from this example optimization formula- 
tion are possible. This method is useful to improve the 
conservatism of the previously described prior-art Q- 
parameterization control design method where the inher- 
ent structure of the uncertainties are commonly ignored 
by simply formulating sigma(H(Q)) which is equivalent 
with a un-scaled H constraint. Fig. 6 shows an example 

00 

flowchart two compare the method in this embodiment 
and the prior-art method. Numerical examples can be 
found in [K. Tsai, Design of Feedforward and Feedback 
Controllers by Signal Processing and Convex Optimization 
Techniques, chapter 3], which is claimed for the invention. 
[0034] |f j n step 5 5 o, the designer chooses to use Q- 

optimization as in step 530, and formulate the robustness 



measure based on a frequency-by-frequency gridding op- 
timization formulation, then in the iteration step 5 40, no 
curve-fitting to the frequency response of D is required; 
and in step 5 20, the frequency response data of the un- 
certainty weights can be provided without curve-fitting it 
to a parametric model. Some prior-art publication ex- 
plaining the procedure to perform this frequency gridding 
optimization can be found in [B. Rafaely et al, "H /H ac- 
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tive control of sound in a headrest: design and implemen- 
tation," IEEE Trans. Control System Technology, vol. 7, no. 
1, Jan. 1999][P. Titterton, "Practical method for con- 
strained-optimization controller design: H or H opti- 

2 oo 

mization with multiple H and/or H constraints," IEEE 
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Proceedings of ASILO 1996]][A. Lanzon et al "A Frequency 
Domain Optimisation Algorithm for Simultaneous Design 
of Performance Weights and Controllers in mu-Synthesis", 
Proceedings of the 38th IEEE Conference on Decision and 
Control, Vol. 5, pp. 4523-4528, Phoenix, AZ, USA, Dec 
1999] [K. Tsai et al, "DQIT: u-synthesis without D-Scale 
Fitting," American Control Conference 2002, pp. 
493-498]. If the robustness measure is formulated based 
on non-frequency gridding approaches, then parametric 
models in step 5 20 and step 5 40 are still required. 



D-Q ITERATION WITH NONPARAMETRIC PLANT UNCERTAINTY WEIGHTS 
INCORPORATED DIRECTLY WITHOUT CURVE-FITTING 

[0035] As another embodiment, when there is no need to incor- 
porated multiple performance objectives with a robust- 
ness measure, but it is desired to reduce the conservatism 
by reducing the modeling error of the plant uncertainty 
from its experimental data, step 450can be substituted 
with "Optimize Q while fixing D, based on frequency grid- 
ding to formulate robustness measures, with or without 
other performance measures". In step 420, the frequency 
response of the plant uncertainty can be provided directly 
from a estimate of the least conservative nonparametric 
model error weight. One of the prior-art publications 
showing a method to perform the estimation is [H. Hindi 
et al, "Identification of Optimal Uncertainty Models from 
Frequency Domain Data," Proceedings of the IEEE Confer- 
ence on Decision and Control, 2002]. It should reduce the 
conservatism due to the use of a simplified and conserva- 
tive model error weighting function of the standard D-K 
iteration. 

D-Q ITERATION WITH ROBUSTNESS MEASURE FORMULATED WITH 
DECISION VARIABLES BEING THE FREQUENCY RESPONSE OF Q 

[0036] when formulating a robustness measure based on fre- 



quency-gridding optimization, it is common practice to 
use the coefficients of the free controller design parame- 
ter Q as the decision variables, and the robustness mea- 
sure objectives or constraints are specified frequency- 
by-frequency with respect to the frequency response of Q. 
The frequency response and the coefficients of the free 
controller design parameter Q can be related by a discrete 
Fourier transform relationship. However, there are cases 
where it is preferred to formulate the optimization prob- 
lem in terms of the frequency response of Q directly. One 
example can be found in [B. Boulet et al An LMI Approach 
to IMC-Based Robust Tunable Control, American Control 
Conference 2003, pp 821]. However, the frequency re- 
sponse of Q can not be optimized independently without 
any constraints, because its inverse discrete Fourier trans- 
form may not be periodically causal, which can cause 
problems when implementing the Q filter. Periodically 
causality is a property of discrete signals which is ex- 
plained in [V. Oppenheim et al, Discrete-Time Signal Pro- 
cessing, Prentice Hall, 1989]. The invention discloses a 
method to impose periodical causality on the frequency 
response of Q. The method is to use Hilbert Transform to 
relate the real part and imaginary part of the frequency 



response of Q at each of a selected set of frequencies. 
One example of the mathematical formulation is found in 
[K. Tsai and H. Hindi, "DQIT: u-synthesis without D-Scale 
Fitting," American Control Conference 2002, pp. 493-498] 
which is claimed for the invention. 
[0037] with this optimization formulation, another embodiment 
to improve the prior-art D-K iteration in Fig 3 is to substi- 
tute step 450with "Optimize Q while fixing D, based on 
frequency gridding to formulate robustness measures, 
with or without other performance measures" and the fre- 
quency gridding formulation is based using the frequency 
response of the free controller design parameter Q as the 
decision variables, with periodical causality constraints 
imposed on the decision variables. 

A COMPARISON OF D-K ITERATION WITH D-Q ITERATION WITHOUT 
D-FITTING 



Table 1 Comparisons of DKIT and DQIT, 6^=271-1 rad/s 



Iteration 


1 


2 


3 


4 


5 


K order 


6 


10 


10 


10 


10 


total D order 


0 


4 


4 


4 


4 


jUmax (DKIT) 


28.10 


18.99 


16.26 


14.93 


14.32 


JUmax (DQIT) 


31.66 


6.78 


3.55 


3.14 


3.13 



[0039] Table 1 shows a numerical example for the advantage of 
using D-Q iteration versus the standard prior-art imple- 
mentation of D-K iteration. In order to make a one- 
to-one comparison, the control design optimization for- 
mulation in step 450 comprises only a robustness mea- 
sure, which the upper bound of u. It shows that both D-K 
iteration and D-Q iteration can reduce the upper bound of 
u monotonically. However, in D-K iteration, to avoid nu- 
merical conditioning problems during the D-scale fitting, 
a lower fitting order has to be used, thus u converges 
slowly at about 14, comparing to the monotonic decrease 
of u down to 3.13 with D-Q iteration without D-scaling 
fitting. The details of this numerical example is described 
in [K. Tsai et al, "DQIT: u-synthesis without D-Scale Fit- 
ting," American Control Conference 2002, pp. 493-498] 
which is claimed for the invention. 



